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ABSTRACT 



The Surface Detector array of the Pierre Auger Observatory can detect neu- 
trinos with energy E v between 10 17 eV and 10 20 eV from point-like sources across 
the sky south of +55° and north of —65° declinations. A search has been per- 
formed for highly inclined extensive air showers produced by the interaction of 
neutrinos of all flavours in the atmosphere (downward-going neutrinos), and by 
the decay of tau leptons originating from tau neutrinos interactions in the Earth's 
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crust (Earth- skimming neutrinos). No candidate neutrinos have been found in 
data up to 2010 May 31. This corresponds to an equivalent exposure of ~3.5 
years of a full surface detector array for the Earth-skimming channel and ~2 years 
for the downward-going channel. An improved upper limit on the diffuse flux of 
tau neutrinos has been derived. Upper limits on the neutrino flux from point- 
like sources have been derived as a function of the source declination. Assuming 
a differential neutrino flux kps ■ E~ 2 from a point-like source, 90% CL. upper 
limits for kps at the level of ~ 5 x 10~ 7 and 2.5 x 10~ 6 GeV cm -2 s -1 have been 
obtained over a broad range of declinations from the searches of Earth- skimming 
and downward-going neutrinos, respectively. 

Subject headings: astroparticle physics — cosmic rays — neutrinos — telescopes 



1. INTRODUCTION 



The nature and production mechanism s of ultra-high energy cosmic rays (UHECRs), 
with e n ergies above 10 18 eV , are st ill unknown iNagano fc Watson I (120001 ) ; iBhattacharjee fc Sigl 



(120001 ); lHalzen fc Hooper I (120021 ) . The observat ion of UHECR s makes an associated flux 



of ultra-high energy cosmic neutrinos (UHEz/s) iBecker I (120081 ) very likely. All models of 
UHECR production predict neutrinos as a result of the decay of charged pions gener- 
ated in interactions of cosmic rays within the sources themselves ( "astrophysical" neutri- 
nos), and/or in their propagation through background radiation fields ( "cosmogenic" neu- 

In fact, charged pions, which are 



( 


1969): 


Stecker 


( 


1973) 



photoproduced by UHECR protons interacting with the Cosmic Microwave Background ra- 
diation, decay into UHEz/s. However, the predicted flux has large uncertainties, since it 
depends on the U HECR spectrum and on the spatial di stribution and cosmological evolu- 



tion of the sources Becker 



(120081 ): lAhlers et al. I (120101 ) 



are heavy nuclei, the UHEz/ yield is strongly suppressed lAve et al. I ( 120051 ). 



Koteraeta 



(|2010h . If UHECRs 



The observation of UHE neutrinos would open a new window to the Universe. Neutri- 
nos travel unaffected by magnetic fields and can give information on astrophysical regions 
that are otherwise hidden from observation by large amounts of matter. The discovery of 
astrophysical neutrino sources would shed light on the long-standing question of the origin 
of cosmic rays, and clarif y the production mechanism of th e GeV-TeV gamma-rays observed 
on Earth iGaisser] Jl995h : klvarez-Muniz fc Halzen I pxjj. 



The Pierre Auger Observatory Abraham et al. I (120041 ) - located in the province of 



Mendoza, Argentina, at a mean altitude of 1400 m above sea level (~875 g cm 



was 
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designed to me asure extensive a i r show ers (EAS) induced by UHECRs. The Fluorescence 
Detector (FD) lAbraham et al. I (j2010al ) comprises a set of imaging telescopes to measure 
the light emitt ed by excited atmosphe ric nitrogen molecules as the EAS develops. A Surface 



Detector (SD) lAUekotte et al. I ( 120081 ). measures EAS particles at ground with an array of 
water-Cherenkov detectors ("stations"). Each SD station contains 12 tonnes of water viewed 
by three 9" photomultipliers (PMTs). Arranged on a triangular grid with 1.5 km spacing, 
1660 SD stations are deployed over an area of ~ 3000 km 2 , overlooked by 27 fluorescence 
telescopes. 

Although the primary goal of the SD is to detect UHECRs, it can also identify ultra-high 
energy neutrinos. Neutrinos of all flavours can interact at any atmospheric depth through 
charged or neutral currents and induce a "downward-going" (DG) shower. In addition, 
tau neutrinos can undergo charged current interactions in the Earth crust and produce a 
tau lepton which, after emerging from the Earth surface and decaying in the atmosphere, 
will induce an "Earth-skimming" (ES) upward-going shower. Even if tau neutrinos are not 
expected to be produced at the astrophysical source, approximately equal fluxes for each 
neutrino flavour should reach the Earth as a result of neutrino oscillations over cosmological 
distances. Neutrino candidate events must be identified against the overwhelming back- 
ground of showers in itiated by standard UHEC Rs (protons or nuclei) and, in a much smaller 
proportion, photons lAbraham et al. I (l2010bl ). Highly inclined (zenith angle 9 > 75°) ES 
and DG neutrino-induced showers will present a significant electromagnetic component at 
the ground ("young" showers), producing signals spread over hundreds of nanoseconds in 
several of the triggered SD stations. Inclined showers initiated by standard UHECRs are, by 
contrast, dominated by muons at ground level ("old" showers), with signals typically spread 
over only tens of nanoseconds. Thanks to the fast sampling (25 ns) of the SD digital elec- 
tronics, several observables sensitive to the signal time structure can be used to discriminate 
between young and old showers, allowing for detection of UHEz/s. Candidates for UHEz/s 
are searched for in inclined showers in the ranges 75° < 9 < 90° and 90° < 9 < 96° for the 
DG and ES analysis, respectively. 



2. LIMITS ON THE DIFFUSE FLUX OF UHE TAU NEUTRINOS 



An upper limit on the diffuse flux of tau neutrinos from the search of Earth-skimming 
ev ents in data through 2 008 April 30 (~2 years of exposure with a full SD array) was reported 



m 



Abraham et al.l (120091 ) . Here, the search is extended to data until 2010 May 31 (~3.5 years 



of exposure with a full SD arr ay), and an improved limit is obtained. A preliminary report 



of this result was presented in iGuardincerri et al. I (120111 ) 
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Details of the neutrino selection procedure, of the calculation o f the detector exposur e 
for ES showers, and of sources of systematic uncertainties are given in lAbraham et al.l ( 120091 ) . 
The neutrino selection criteria were optimized with an early data set collected between 2004 
November 1 and 2004 December 31. By using data rather than Monte Carlo simulations, 
all possible detector effects and shower-to-shower fluctuations, which constitute the main 
background to UHEz/s and may not be well reproduced by simulations, are taken into ac- 
count. The neutrino selection established with the training sample was then applied to a 
"blind search sample" of data collected between 2004 January 1 and 2010 May 31 (exclud- 
ing November and December 2004). The blind search sample is equivalent to ~ 3.5 years 
of data collected by a fully efficient SD array, i.e. with all stations working continuously. 
The time evolution of the SD array, which was growing during the construction phase, as 
well as the dead times of individual stations, were accounted for in this calculation. The 
integrated exposure for detection of ES tau neutrinos as a function of energy is shown in 
Figure [H No neutrino candidates were found in the blind search. Assuming a differential 
spectrum <&{E V ) = dN v /dE v = k ■ E~ 2 for the diffuse flux of UHEz/s and zero background 



Abraham et al.l (120091 ); lAbreu et al.l (j2011al ). a 90% C.L. upper limit on the integrated flux 



of tau neutrinos is derived: 



k < 3.2 x 10 



-8 



GeV cm 



2 s- 1 sr- 1 . 



Systematic uncertainties in the exposure were taken into account in the upper li mit by using a 

semi-B ayesian extension IConrad et al. I (120031 ) of the Feldman- Cousins approach iFeldman &: Cousins 
(119981 ). The limit, shown as a horizontal line in Figure [2J is valid in the energy range 
1.6 x 10 17 eV < E v < 2.0 x 10 19 eV, where ~ 90% of neutrino events would be detected for 
a E~ 2 flux . Also shown is the 90% C.L. upper limit in differential form, where the limit 
was calculated independently in each energy bin of width 0.5 in log ^ E v . The in t egrated 
and differential limits from the search for downward-going neutrinos lAbreu et al.l ( l2011af ) 
at the Pierre Auger Observatory, based on a "blind search sample" of data collected from 
2007 November 1 until 2010 May 31 (equivalent to ~ 2.0 years of exposure with the full 
SD array) , are also shown in Fig ure |2j together with limits fro m the IceCube Neutri no Ob- 
servatory lAbbasi et al. I (j2011al ) and the ANITA experiment iGorham et al. I (120101 ) . The 
shaded area in Figure [2] brackets the cosmogenic neutrinos fluxes predicted under a wide 
range of assumptions for the cosmological evolution of the sources, for the transition be- 
tween th e galactic and extrag alactic component of cosmic rays, and for the UHECR com- 
position iKotera et al. I ( 120101 ). The corresponding number of cosmogenic neutrino events 
expected in the blind search sample ranges bet ween 0.1 and . 3, ap proximately. For the 
diffuse flux of cosmogenic neutrinos predicted in lAhlers et al. I (120101 ) . 0.6 neutrino events 
are expected at the Pierre Auger Observatory with the integrated exposure of the present 
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analysis, to be compared wit h 0.43 events expected in the 333.5 days of live-time of the 



IceCube-40 neutrino telescope lAbbasi et al. J (2011a) . Th e curre nt bound to a cosmogenic 



neutrino flux with energy dependence as in lAhlers et al. I (120101 ) and shown in Figure [2] is 
4 times larger than the predicted value. With the current selection criteria the exposure 
accumulated in ~10 more years with the Pierre Auger Observatory may exclude this cosmo- 
genic neutrino flux at 90% C.L.. Notice that the maximum sensitivity of the Pierre Auger 
Observatory, obtained for E u ~ 10 18 eV, matches well the peak of the expected neutrino 
flux. 



3. SENSITIVITY TO POINT-LIKE SOURCES 

The neutrino search at the Pierre Auger Observatory is limited to highly inclined show- 
ers, with zenith angles between 90° and 96° in the Earth-skimming analysis, and between 
75° and 90° in the downward-going analysis. Thus, at each instant, neutrinos can be de- 
tected only from a specific portion of the sky corresponding to these zenith angle ranges. A 
point-like source of declination 5 and right ascension a (equatorial coordinates) is seen at 
our latitude (A = —35.2°), at a given sidereal time t, with a zenith angle 9(t) given by: 

cos 6(t) = sin A sin 5 + cos A cos 5 sin(27rf/T — a) , (2) 

where T is the duration of one sidereal day. From equation [2J the fraction of a sidereal 
day during which a source is detectable (i.e., within the zenith angle ranges for the ES and 
DG analyses) is shown in Figure [3j it depends only on the source declination. The SD of 
the Pierre Auger Observatory is sensitive to point-like sources of neutrinos over a broad 
declination range spanning north of S ~ —65° and south of 5 ~ 55°. The regions of the sky 
close to the Northern (5 = 90°) and Southern (5 = —90°) Terrestrial Poles are not accessible 
by this analysis. As an example, Centaurus A (5 ~ —43°) is observed ~ 7% (~ 15%) of one 
sidereal day in the range of zenith angles corresponding to the ES (DG) search. The peaks 
in Figure [3] are a consequence of the relatively smaller rate of variation of zenith angle with 
time for directions near the edges of the range accessible to this analysis. 

The exposure of the SD as a function of the neutrino energy and of the source position 
in the sky, £(E U ,5), is evaluated by folding the SD aperture with the neutrino interaction 
probability and the selection efficiency for each neutrino channel. The proce dure is identical 



to tha t used for the c alculat io n of the exposure for a diff use flux of UHEz/s ([Abraham et al. 



(2009 


); 


Abreu et al. 


( 


2011a 


); 



angle integration over the sky. The integration over the blind search time period takes 
into account the growth of the SD array during its construction phase and dead times of 
individual stations. For example, the exposure for the DG analysis is given by: 
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1 E 




cos0(t) 8i(f,E u ,e(t),D,t) dA dD dt 



(3) 



where the integration is performed over the area A of the SD, the interaction depth D 
of the neutrino, and the search peri od. In equation [31 m is the ma ss of a nucleon, ai(E v ) 
is the neutrino-nucleon cross-section ICooper-Sarkar fc Sarkar I (120081 ). and Si is the neutrino 
selection efficiency, with the sum running over the 3 neutrino flavours (a;, = 1, corresponding 
to a 1:1:1 flavour ratio) and over the neutrino charged and neutral current interactions. The 
dependence of e on several parameters (the point of impact at ground of the shower core, 
f, the neutrino interaction depth, its energy and zenith angle, and time) is also explicitly 
included in equation [3j The dependence of the exposure on the source declination comes 
through 8(t) as obtained from equation [2j When integrating over time, only those periods 
when the source is within the zenith angle range of the neutrino selection are considered. 
The exposure for ES neutrinos is derived analogously to equation [3J 

Changes in the detector configuration during data taking, due to the dead times of 
the SD stations, and to the increase of the array size during the construction phase, may 
introduce a dependence of the exposure on the right ascension. In particular, fluctuations 
in the number of stations cause a small diurnal variation, but this effect is only apparent 
in solar time. When averaged over a large number of sidereal days, as in this ana l ysis, th e 
modulation in right ascension caused by this effect is less than 1% lAbreu et al. I ( 1201 lbl ). 
For this reason, the dependence of the exposure on a has been neglected in the evaluation 
of the upper limits. 

Due to the finite resolution of the SD on the reconstruction of the variables used in the 
selection of neutrino-induced showers, events close to the edges of the zenith angle range for 
the neutrino selection may be wrongly rejected (or wrongly accepted). In the exposure as 
given in equation [3] we account for this effect by evaluating the selection efficiency e through 
Monte Carlo simulations. 

Several other sour c es of systematic unce rtainties on the exposure have been investigated 
( 1 Abraham et al.l ( 120091 ): lAbreu et al.l fcOllaf )). For the DG analysis, the major contributions 
in terms of deviation from a reference exposure come from the knowledge of neutrino-induced 
shower simulations (+9%, —33%), of the neutrino cross-section (±7%), and of the topography 
(±6%). Only uncertainties c ompatible with the conventional NLO DGLAP formalism of v 
cross-section calculation - see lCooper-Sarkar Sarkar I ( 120081 ) for details - have been consid- 
ered. We have not accounted for gluon saturation models t hat would give rise to considerabl e 
smaller v cross- sect ions (as small as a factor ~ 2 at 10 18 eV lHenley &: Jalilian-Marian I ( 120061 ): 
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Armesto et al.l ( 120081 )). and hence to a larger systematic uncertainty than the one quoted 
here. For the ES analysis, the systematic uncertainties are dominated by the tau energy losses 
(+25%, -10%), the shower simulations (+20%, -5%) and the topography (+18%, 0%). 



4. LIMITS ON THE FLUX OF UHE NEUTRINOS FROM POINT-LIKE 

SOURCES 

The expected number of neutrino events in an energy range [E m i n , -^max] from a point-like 
source located at a declination 5 is given by: 

/■E max 

^ e P x°p" e W = / F(K) S(E U , 5) dE u , (4) 

where F(E V ) is the flux of UHEz/s from the source. No candidate events were selected using 
the .E'S' and DG analyses. Under the conservative assumption of zero background, a 90% C.L. 
upper limit on the neutrino flux from point-like sources is derived. To set the upper limit, a 
differential flux F(E U ) = kps(5) ■ E~ 2 was assumed, as well as a 1:1:1 neutrino flavour ratio. 
Systematic uncertainties on the exposure were calculated using the semi-Bayesian approach 
described above in section 2. 

In Figure HI the 90% C.L. upper limits on kps derived from the ES and DG analyses 
are shown as a function of source declination. Limits for kps at the level of ~ 5 x 10~ 7 and 
2.5 x 10~ 6 GeV cm" 2 s" 1 were obtained over a broad range of declinations from the searches 
of Earth-skimming and downward-going neutrinos, respectively. 

The shape of the declination-dependent upper limits is largely determined by the fraction 
of time a source is within the field of view of the ES or DG analyses (cf. Figure [3]), and, to 
a lesser extent, by the zenith angle dependence of the exposure. 

The upper limits are derived for neutrinos in the energy range 1.6 x 10 17 eV— 2.0 x 10 19 eV 
for the Earth-skimming analysis, and in the energy range 1 x 10 17 eV — 1 x 10 20 eV for the 
downward-going analysis, with a negligible dependence of these energy intervals on the source 
declination. These are the best limits around 1 EeV. 

The IceCube Neutrino Observatory and the Antares Neu trino Telescope have also 



search ed for UHEz/s from point-like sources (lAbbasi et al. I (1201 Ibt ) and lAdrian-Martmez et al 



(120 111 ), respectively). The bounds obtained by these two experiments apply to energies below 



the Auger energy range. 



Limits for the particular case of the active galaxy Centaurus A, a potential source of 
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UHECRs, are shown in Figure O together with constraints from other e xperiments. The pre- 



dicted fluxes for two theoretical models of UH Ez/ production -in the jetslCuoco fc Hannestad 



( 120081 ) and close to the core of Centaurus A iKachelriess et al. I ( 120091 ) - are also shown for 
co mparison. The expec t ed nu mber of events in our blind search samples for a flux like 



m 



Cuoco fc Hannestad I (120081 ) is about 0.1 and 0.0 2 for t he ES and DG selection respec- 

( 120091 ) being one order of magnitude 



tively, the expected number for IKachelriess et al. 
smaller. 



5. SUMMARY 

The sensitivity of the Pierre Auger Observatory to point-like sources of neutrinos with 
ultra high-energy has been obtained. Highly inclined extensive air showers produced by the 
interaction of neutrinos of all flavours in the atmosphere and by the decay of tau leptons 
originating from tau neutrinos interactions in the Earth crust were searched for, and discrim- 
inated from the background of standard UHECRs thanks to the distinctive time structure of 
the signals measured by the Surface Detector array. The search for neutrinos was performed 
over a broad range of declination, north of ~ —65° and south of ~ 55°, and for neutrino 
energies between 10 17 eV and 10 20 eV. 

No neutrino candidates were found in data collected through 2010 May 31, and an im- 
proved upper limit on the diffuse flux of tau neutrinos was correspondingly placed. Also, the 
exposure for neutrino fluxes from point-like sources was evaluated as a function of source 
declination. Upper limits at 90% CL. for neutrino fluxes from point-like sources were es- 
tablished, which are currently the most stringent at energies around and above 1 EeV in a 
large fraction of the sky spanning more than 100° in declination. 
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Fig. 1. — Exposure of the Surface Detector of the Pierre Auger Observatory for Earth- 
skimming neutrino initiated showers as a function of the neutrino energy, for data collected 
between 2004 January 1 and 2010 May 31. 
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Single flavour neutrino limits (90% CL) 
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Fig. 2. — Differential and integrated upper limits at 90% CL. on the single ffavour E~ 2 
neutrino flux from the search for downward-going and Earth-skimming neutrinos at the 
Pierre Auger Observatory. Integrated upper limits are indicated by horizontal lines, with 
the corresponding differential limits being repr esented by segments of width 0.5 in \og w E v . 
Limits from the IceCube Neutrino Observatory lAbbasi et al. I (j2011al ) and from the ANITA 
experiment iGorham et al. I ( 120101 ) are also shown after proper rescaling to account for sin- 
gle flavour neutrino flux and different ener g y bin ni ng. Prediction s for c osmogenic neutrinos 
under different assumptions I Ahlers et al. I ( 120101 ) ; iKotera et al. I (120101 ) are also shown, al- 
though predictions almost one order of magnitude lower or higher exist. 



-17- 




Source declination 8 [deg] 

Fig. 3. — Fraction of a sidereal day having a point-like source at declination 5 detectable 
by the Pierre Auger Observatory with the Earth-skimming and downward-going neutrino 
selection. 
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Fig. 4. — Upper limits at 90% C.L. on a single flavour E~ 2 flux from a specific point-like 
source as a function of the source declination. The bounds from the Earth-skimming and 
downward-going neutrino analyses hold for a neutrino energy range 10 17 — 10 20 eV (see text 
for details). 
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Centaurus A - Single flavour neutrino limits (90% CL) 
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Fig. 5. — Upper limits at 90% CL. on a single flavour E~ 2 flux from the active galaxy 
Centaurus A from the Earth-skimming and downward- going neutrino analyse s, together 
with bo unds from the IceC ube Neutrino Observatory lAbbasi et al. I (I2011bl ) and LU- 
NASKA^ 



James et al. 



(120111) . T he predictions for two models of UHEz/ production - in the 



j ets ICuoco fc Hannestad I (120081 ). and close to the core of Centaurus A Kachelriess et al. 
J200J - are also shown. 



-20 - 



REFERENCES 

Abbasi, R., Abdou, Y., Abu-Zayyad, T., et al. (IceCube Collaboration) 2011a, Phys. Rev. 
D, 83, 092003 

Abbasi, R., Abdou, Y., Abu-Zayyad, T., et al. (IceCube Collaboration) 2011b, ApJ, 732, 18 

Abraham, J., Abreu, P., Aglietta, M., et al. (Pierre Auger Collaboration) 2009, Phys. Rev. 
D, 79, 102001 

Abraham, J., Abreu, P., Aglietta, M., et al. (Pierre Auger Collaboration) 2010a, Nucl. 
Instrum. Methods Phys. Res. A, 620, 227 

Abraham, J., Abreu, P., Aglietta, M., et al. (Pierre Auger Collaboration) 2010b, Astropart. 
Phys., 31, 399 

Abraham, J., Aglietta, M., Aguirre, I. C, et al. (Pierre Auger Collaboration) 2004, Nucl. 
Instrum. Methods Phys. Res. A, 523, 50 

Abreu, P., Aglietta, M., Ahn, E. J., et al. (Pierre Auger Collaboration) 2011a, Phys. Rev. 
D, 84, 122005 

Abreu, P., Aglietta, M., Ahn, E. J., et al. (Pierre Auger Collaboration) 2011b, Astropart. 
Phys., 34, 627 

Adrian-Martinez, S., Aguilar, J. A., Samarai, I. A., et al. (Antares Collaboration) 2011, ApJ, 
743, L14 

Ahlers, M., Anchordoqui, L.A., Gonzalez- Garcia, M.C., Halzen, F., k Sarkar, S. 2010, As- 
tropart. Phys., 34, 106 

Allekotte, I., Barbosa, A. F., Bauleo, P., et al. 2008, Nucl. Instrum. Methods Phys. Res. A, 
586, 409 

Alvarez-Muhiz, J., k Halzen, F. 2002, ApJ, 576, 33 

Armesto, N., Merino, C, Parente, G., k Zas, E. 2008, Phys. Rev. D, 77, 013001 

Ave, M., Busca, N., Olinto, A. V., Watson, A. A., k Yamamoto, T. 2005, Astropart. Phys., 
23, 19 

Becker, J.K. 2008, Phys. Rep., 458, 173 

Berezinsky, V.S., k Zatsepin, G.T. 1969, Phys. Lett. B, 28, 423 



- 21 - 



Bhattacharjee, P., k Sigl, G. 2000, Phys. Rep., 327, 109 

Cooper-Sarkar, A., k Sarkar, S. 2008, J. High Energy Phys., JHEP01(2008)075 

Conrad, J., Botner, O., Hallgren, A., k Perez de los Heros, C. 2003, Phys. Rev. D, 67, 012002 

Cuoco, A., k Hannestad, S. 2008, Phys. Rev. D, 78, 023007 

Feldman, G.J., k Cousins, R.D. 1998, Phys. Rev. D, 57, 3873 

Gaisser, T.K., Halzen, F., k Stanev, T. 1995, Phys. Rep., 258, 173 

Gorham, P., Allison, P., Baughman, B. M., et al. (ANITA Collaboration) 2010, Phys. Rev. 
D, 82, 022004 (erratum: Phys. Rev. D, 85, 049901(E)) 

Guardincerri, Y. et al. (Pierre Auger Collaboration) 2011, in Proc. 32 st International Cosmic 
Ray Conference 2011 (Contribution No. 0682), Beijing, China 

Halzen, F., k Hooper, D. 2002, Rep. Prog. Phys., 65, 1025 

Henley, E. M., k Jalilian-Marian, J. 2006, Phys. Rev. D, 73, 094004 

James, C.W., Protheroe, R. J., Ekers, R. D., et al. (LUNASKA Collaboration) 2011, MN- 
RAS, 410, 885 

Kachelriess, M., Ostapchenko, S., k Tomas, R. 2009, New J. Phys., 11, 065017 
Kotera, K., Allard, D., k Olinto, A. 2010, J. Cosmol. Astropart. Phys., JCAP10(2010)013 
Nagano, M., k Watson, A. 2000, Rev. Mod. Phys., 72, 689 
Stecker, F.W. 1973, Astrophys. Space Sci., 20, 47 



This preprint was prepared with the A AS IATgX macros v5.2. 



